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Abstract 

NGC 6388 and NGC 6441 are two massive Galactic bulge globular clusters which share many properties, including 
the presence of an extended horizontal branch (HB), quite unexpected because of their high metal content. In this paper 
we use HST's WFPC2, ACS, and WFC3 images and present a broad multicolor study of their stellar content, covering 
all main evolutionary branches. The color-magnitude diagrams (CMDs) give compelling evidence that both clusters host 
at least two stellar populations, which manifest themselves in different ways. NGC 6388 has a broadened main sequence 
(MS), a split sub-giant branch (SGB), and a split red giant branch (RGB) that becomes evident above the HB in our data 
set; its red HB is also split into two branches. NGC 6441 has a split MS, but only an indication of two SGB populations, 
while the RGB clearly splits in two from the SGB level upward, and no red HB structure. The multicolor analysis of 
the CMDs confirms that the He difference between the two main stellar populations in the two clusters must be similar. 
This is observationally supported by the HB morphology, but also confirmed by the color distribution of the stars in the 
MS optical band CMDs. However, a MS split becomes evident in NGC 6441 using UV colors, but not in NGC 6388, 
indicating that the chemical patterns of the different populations are different in the two clusters, with C, N, O abundance 
differences likely playing a major role. We also analyze the radial distribution of the two populations. 



Subject headings: 
motions 

1. Introduction 



globular clusters: individual (NGC 6388, NGC 6441) - photometry - Stars: Population II - C-M diagrams - proper 



Since the 1990s, NGC 6388 (/ = 345?56, b = -6?74) 
and NGC 6441 (/ = 353?53, b = -5?01) have been consid- 
ered as "twin" clusters because of many common features. 
They are both located in the bulge of our Galaxy, they are 
two of the most luminous and massive (~1.6xl0'' M©, Pryor 
& Meylan 1993) globular clusters (GCs) of the Milky Way, 
and they are both metal rich ([Fe/H]=-0.55 for NGC 6388 and 
[Fe/H]=-0.46 for NGC 6441, Harris 1996, 2010 edition). De- 
spite their high metal content, their most intriguing common 
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property is an extended blue horizontal branch (HB, Rich et al. 
1997), which makes them among the most extreme examples 
of the so-called second-parameter problem. Moreover, the red 
clump (RC) of the HB is noticeably tilted, with bluer stars be- 
ing brighter (in the optical bands) than redder ones (Raimondo 
et al. 2002). 

Since the work of Rich et al. (1997), attempts to understand 
what is special about these two clusters have focused on these 
two peculiarities. Sweigart & Catelan (1998) first proposed 
high helium as the origin of the blue HB extension, in a sce- 
nario in which the envelope of fast-rotating stars, while ascend- 
ing the red giant branch (RGB), would be enriched in helium 
by meridional currents. Others, focusing on the HB tilt, argued 
for differential reddening (DR) as the cause of the tilt itself 
(e.g., Raimondo et al. 2002), but even though these clusters are 
certainly affected by DR, the size of the tilt is too large to be 
explained by DR alone, and finally, DR has been excluded by 
Busso et al. (2007), who found that the HB tilt persisted in all 
parts of their field, in each of the two clusters. 

It became evident that traditional evolutionary scenarios 
could not account for the properties of these two clusters fol- 
lowing the discovery that their RR Lyrae variables have ex- 
ceptionally long periods (something that could be explained 
by high helium but certainly not by DR). With an average 
fundamental-mode period of (Pab) = 0?71 for NGC 6388, and 
0?76 for NGC 6441, these two clusters were legitimately con- 
sidered to be the prototype of a third Oosterhoff type (the av- 
erage periods of Oosterhoff I and Oosterhoff II clusters being 
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Table 1 

Log of HST Observations for NGC 6388 



GO 


PI 


Instrument 


Filter 


Exposures 


Ar^ 


Epoch 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


9835 


G. Drukier 


ACS/HRC 


F555W 


51xl55s 


o;o5 


30 Oct 2003 








F814W 


5x25 s, 2x469 s, 10x505 s 


o:o5 


30 Oct 2003 


10350 


H. Cohn 


ACS/HRC 


F330W 


2xl266s4xl314s 


o:o5 


7 Apr 2006 








F555W 


3x155 s 


o:o5 


7 Apr 2006 






ACS/WFC 


F660N 


2x1123 s, 4xll71s 


4;05 


7 Apr 2006 


10474 


G. Drukier 


ACS/HRC 


F814W 


5x25 s, 8x501 s, 4x508 s 


o:o5 


4-8 Apr 2006 


9821" 


B. Protz 


ACS/WFC 


F435W 


6xlls 


0:08 


2003-2004 








F555W 


6x7 s 


0:08 


2003-2004 


10775 


A. Sarajedini 


ACS/WFC 


F606W 


Ix40s,5x340s 


o;27 


6 Apr 2006 








F814W 


1x40 s, 5x350 s 


0:27 


7 Apr 2006 


11739 


G. Piotto 


ACS/WFC 


F475W 


4x788 s 


3:86 


3 Jul 2010 








F814W 


1x67 s, 4x155 s 


4:02 


30Jun 2010 






WFC3/UVIS 


F390W 


6x880 s 


2:35 


3 Jul 2010 






WFC3/IR 


F160W 


6x199 s, 4x249 s 


2;i8 


30 Jun 2010 



t Average distance of data set from cluster center. 
"SNAP program. 



Table 2 

Log of HST Observations for NGC 6441 



GO 


PI 


Instrument 


Filter 


Exposures 


Ar"f 


Epoch 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


5667 


B. H. Margon 


WFPC2 


F336W 


2x500 s 


0:05 


8 Aug 1994 


8251 


H. A. Smith 


WFPC2 


F439W 


3x50 s, 15x80 s 


0:10 


1999-2000 








F555W 


18xlOs 


0:01 


1999-2000 


8718 


G. Piotto 


WFPC2 


F336W 


lx260s, 2x400 s 


0:02 


16 Sep 2001 


9835 


G. Drukier 


ACS/HRC 


F555W 


36x240 s 


o:o7 


1 Sep 2003 








F814W 


5x40 s, 1x413 s, 10x440 s 


0:07 


1 Sep 2003 


10775 


A. Sarajedini 


ACS/WFC 


F606W 


1x45 s, 5x340 s 


0:27 


25 May 2006 








F814W 


1x45 s, 5x350 s 


0:27 


25 May 2006 


11739 


G. Piotto 


ACS/WFC 


F475W 


3x768 s, 3x803 s 


3:86 


4-8 May 2010 










1x765 s, 1x800 s 


3:86 


30 May 2011 








F814W 


2x67 s, 6x155 s 


4:02 


4-8 May 2010 










2x155 s 


4:02 


30 May 2011 






WFC3/UVIS 


F390W 


9x884 s 


2:35 


4-8 May 2010 










3x884 s 


2:35 


30 May 2011 






WFC3/IR 


F160W 


8x199 s, 8x249 s 


2:18 


7 May 2010 



t Average distance of data set from cluster center. 



0?55 and 0?65, respectively, Layden et al. 1999; Pritzl et al. 
2000, 2001,2002, 2003). 

The discovery of multiple stellar populations with enhanced 
helium in several GCs (see Piotto 2009 for a recent review) 
offers a new perspective for understanding the many apparent 
anomalies that these two clusters share. After it was shown that 
high helium could be the cause of both the blue HB extension 
and the unusual period of the RR Lyrae (D'Antona & Caloi 
2004; Caloi & D'Antona 2007; Busso et al. 2007), the question 
became whether high helium was the result of deep-mixing in a 
minority of RGB stars (as originally suggested by Sweigart & 
Catelan 1998), or was due to the presence of a second, helium- 
enriched stellar population. In the former case no splitting of 



the main sequence (MS) was expected, whereas a sizable split 
should exist if the blue HB and RR Lyrae stars belong to a 
distinct, helium-rich sub-population. 

Evidence for different stellar generations has already been 
found in both clusters with the detection of light-element star- 
to-star variations and the Na-O anticorrelation (Gratton et al. 
2007; Carretta et al. 2007, 2009a). Moreover, NGC 6388 also 
exhibits a split sub-giant branch (SGB, Piotto 2008; Moretti et 
al. 2009; Piotto et al. 2012). Thus, if the blue HB and RR Lyrae 
stars belong to helium-enriched populations in these clusters, 
we should be able to see the effect of enhanced He on the MSs 
(a spread or a split), to quantify the enhancement, and to check 
whether it is consistent with the amount of helium enhance- 
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Fig. 1. — Finding charts for NGC 6388 (a) and NGC 6441 (b) superimposed on star maps from the 7-band 2MASS catalog. Data 
sets are color-coded according to the GO program. There is no ambiguity in using the same color for different cameras (e.g., for 
GO-11739 WFC3/UVIS, WFC3/1R, and ACS/WFC) because of the different size of each camera's FoV (which are shown, to scale, on 
the right edge of the figure). ACSAVFC is the larger 2-CCD footprint, and WFC3/UVIS is the smaller one; the larger 1-CCD footprint 
is WFC3/1R, and the smaller one is ACS/HRC (always placed at the cluster center). Finally, WFPC2 has its own distinctive L-shaped 
footprint, and is used here only for NGC 6441. We give both right ascension and declination (lower abscissa and left ordinate) and 
components relative to the cluster center in arcmin (upper abscissa and right ordinate). The inner dotted circle marks the half-light 
radius r^, while the outer one marks the tidal radius (rt, from Harris 1996, 2010 edition). The core radius (not marked) falls within the 
ACS/HRC pointings. 



ment required to account for the HB properties. 

To this end, we submitted a Hubble Space Telescope (HST) 
proposal (GO-11739, PI G. Piotto) to perform accurate pho- 
tometry of NGC 6388 and NGC 6441 with the Wide Field 
Camera 3 (WFC3). In the present paper, we combine GO- 
11739 data with a large set of archival images from the Ad- 
vanced Camera for Surveys (ACS) and the Wide-Field Plane- 
tary Camera 2 (WFPC2), in order to thoroughly map the mul- 
tiple stellar populations of these two clusters from their MS all 
the way to the SGB, RGB and HB. 

This paper is organized as follows: in Section 2 the vast 
set of HST data used for the project is presented, along with 
the reduction procedures. In Section 3 distinct stellar popula- 
tions are identified in all the color-magnitude-diagram (CMD) 
sequences of the two clusters, and in Section 4 their radial dis- 
tributions within the clusters are presented. In Section 5 we 
estimate the the amount of helium enhancement AF required 
to account for the size of the MS split, and discuss its implica- 
tions at the light of our new observational findings. 

2. Data sets and reductions 

This work is based on proprietary and on archival HST im- 
ages. We made use of exposures taken with three diff'erent 
cameras: (1) ACS, with both the Wide-Field Channel and the 



High-Resolution Channel (WFC and HRC); (2) WFC3, with 
both the Ultraviolet- Visible (UVIS) and the Infrared (IR) chan- 
nels; and (3) WFPC2. A list of the observations is given in 
Table 1 for NGC 6388 and in Table 2 for NGC 6441. 

To summarize briefly: ACS/WFC observations come from 
programs SNAP-9821, GO-10775 andGO-1 1739forNGC 6388; 
from GO-10775 and GO-11739 for NGC 6441. ACS/HRC 
observations are from programs GO-9835, GO-10350, and 
GO- 10474. WFPC2 exposures come from GO-5667, GO- 
8251, and GO-8718. Much of our material is archival, but 
as an essential part of this study we specifically designed the 
new program GO-11739 (PI: Piotto) to detect multiple stel- 
lar populations along the MSs of NGC 6388 and NGC 6441 
by exploiting the power of the new WFC3 cameras. We used 
both WFC3/UV1S and WFC3/1R, in each case taking parallel 
exposures with ACS/WFC. 

Figure 1 shows the footprints of all exposures that are listed 
in Tables 1 and 2, superimposed on stars from the 2MASS J- 
band catalog. A distinctive color is used for each program and 
a distinctive shape and size for each camera. 

Figure 1 and Tables 1 and 2 are crucial for this paper, as 
they exhibit the overlaps of the different data sets, which gov- 
ern much of what follows. Careful study of the figure, along 
with the tables, will enable the reader to see possible filter com- 
binations that can be used to build CMDs at different radial 
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locations in the clusters. Moreover, since some overlapping 
data sets encompass different epochs, one can see where proper 
motions are measurable, information that will prove invaluable 
for selection of cluster members. The reader should particu- 
larly note, for each cluster, the overlap between GO-11739 and 
GO-10775 (much smaller for NGC 6388 than for NGC 6441). 

In this paper we will focus our attention on the more cen- 
tral fields, leaving the analysis of the outer fields to a later pa- 
per This choice will not affect our results, since the outer field 
has been imaged only with filters not designed to highlight the 
multi-population phenomenon along the MS of these two clus- 
ters. Moreover, the outer field does not have enough RGB and 
SGB stars for a study of their radial distribution. 

2.1. ACS/WFC and HRC 

All ACS/WFC _£lt images (i.e., images that were dark- 
and bias-subtracted and flat-fielded, but not resampled) were 
corrected for charge-transfer deficiencies by using the proce- 
dure and the software developed by Anderson & Bedin (2010). 

Photometry and astrometry were carried out with the soft- 
ware tools described by Anderson et al. (2008). Briefly, all ex- 
posures belonging to a specific program were analyzed simul- 
taneously to generate an astrophotometric catalog of stars in 
the field of view (FoV). Stars were measured independently in 
each image by using the spatially varied point-spread function 
(PSF) library models from Anderson & King (2006b), along 
with a spatially constant perturbation specifically derived for 
each exposure, to compensate for small focus changes due to 
spacecraft "breathing". Star positions were corrected for geo- 
metric distortion using Anderson & King (2006b) and Ander- 
son (2007). 

Photometry was calibrated into the Vega-mag flight system 
following recipes in Bedin et al. (2005), and using the zero 
points given in Sirianni et al. (2005). 

The measurement of stellar fluxes and positions in each 
ACS/HRC image was performed by using the publicly avail- 
able routine img2xym_HRC, library PSFs, and the distortion 
correction described in Anderson & King (2006a). We cor- 
rected the photometric zero points in the same way as for the 
WFC photometry. 

2.2. WFPC2 

The WFPC2 images are not _£lt-type (introduced only for 
later instruments), but rather the earlier equivalent _c0f . We 
reduced them with the algorithms described in Anderson & 
King (2000). The field was calibrated to the photometric Vega- 
mag flight system of WFPC2 according to the prescriptions of 
Holtzmanet al. (1995). 

2.3. WFC3/UVIS and IR 

For the WFC3 exposures we used the _flt images. Star 
positions and fluxes on WFC3/UVIS images were measured 
with software based mostly on img2xym_WFI (Anderson et al. 
2006), which implements a set of spatially variable empirical 
PSFs for each individual exposure. We corrected positions for 
geometric distortion using the recipes given in Bellini & Bedin 



(2009) and Beflini, Anderson & Bedin (201 1). Photometry was 
calibrated as in Bedin et al. (2005). 

Infrared images were reduced with software based on that 
of Anderson & King (2006a) for ACS/HRC. We constructed an 
empirical F160W "effective PSF" using the 3-step procedures 
devised by Anderson & King (2000). A detailed description of 
the software, PSF derivation, and geometric distortion solution 
will be presented elsewhere. Photometry was calibrated as in 
Bedin et al. (2005). 

2.4. Sample selections 

Since our aim is fine photometric discrimination, it was a 
great advantage to select the best-measured stars in our sam- 
ples. To accomplish this, we closely followed the recipes and 
used the same selection criteria as given in Milone et al. (2009), 
which makes use of several diagnostics: photometric and as- 
trometric rms, PSF-fit residuals, and amount of scattered light 
from neighboring stars. (This last diagnostic is available only 
for ACS/WFC catalogs because, as we stated above, different 
codes where used for different cameras, and only the one for 
the ACS/WFC provides the additional diagnostic.) 

For more information about our selection procedures we re- 
fer the reader to Milone et al. (2009). 

2.5. Differential reddening 

Both NGC 6388 and NGC 6441 are in the direction of the 
Galactic bulge, where extinction is a serious problem {E{B-V) 
0.37 and 0.47 respectively, according to Harris 1996, 2010 edi- 
tion). Such large reddenings can be expected to be accompa- 
nied by serious differential reddening, which we treated with 
particular care. 

As a first step, for each CMD we computed the direction of 
the reddening vector For this we used the extinction coeffi- 
cients in Tables 14 and 15 of Sirianni et al. (2005) for G2 stars. 
For WFC3 F390W, which is not in those tables, we interpo- 
lated lineady between ACS/HRC F330W and F435W, getting 
AF39ow/£'(fi - V) = 4.49. For the F160W filter we adopted the 
value AFi6nw/£'(B - V) ^ 0.57 (from YES, the York Extinc- 
tion Solveo, and used the Fitzpatrick (1999) extinction law. 
As we did for the WFPC2, we took the extinction values from 
Schlegel, Finkbeiner & Davis (1998). 

We dealt with differential reddening (DR) by means of a 
method that is explained in detail in Section 3.1 of Milone et 
al. (2012a), to which we refer the reader for details; here we 
only describe the method briefly. 

The overall method consists of iterating a core procedure 
which itself has five steps: (1) We rigidly rotate the CMD to 
make the reddening direction parallel to the x axis, so that DR 
becomes a simple shift in x, and the x distance of a star from 
the now- vertical fiducial becomes an estimate of its DR value. 
(2) We next choose a set of "reference stars" — stars whose 
DR values are best suited to serve as standards. Ideally these 
stars should lie in the parts of the CMD where the sequences 



http ://www3 . cadc-ccda.hia-iha.nrc-cnrc.gc . ca/community/ 
YorkExtinctionSolver/ 
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NO CORRECTION 



FIRST ITER 



EINAL CORRECTION 




Fig. 2. — Illustrating the correction for differential reddening. The first panel shows the uncorrected /izpguw vs. mpsgow - '«f814w CMD 
of NGC 6441 around the SGB region. Stars used as a reference for the first-guess DR correction are marked in red, and the fiducial 
sequence in green. The arrow shows the direction of the reddening vector. After a trial correction (second panel), the sequences in the 
CMD are significantly narrower, and we are able to refine our choice of reference stars and adjust the fiducial sequence a little. In the 
final panel is the result of three iterations. Note the sharpness of the SGB, and the splits in the MS and RGB. 



ALL STARS MEMBERSHIP GC MEMBERS FIELD STARS 




Fig. 3. — Illustrating the proper-motion separation for NGC 6388. (a) CMD for stars in common between GO-10775 and GO-11739. 
(b) Sizes of the proper-motion displacements (in pixel/yr); the red line (drawn by hand) shows our membership criterion, (c) CMD of 
proper-motion-selected cluster members; (d) stars rejected as field stars (even though this means the loss of a few members). 
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ALL STARS MEMBERSHIP GC MEMBERS FIELD STARS 




Fig. 4. — Same as Fig. [3] but for NGC 6441. Although the field contamination is much stronger, we have nevertheless selected an 
almost clean sample of cluster members. 



are almost perpendicular to the reddening vector, i.e., the SGB 
and the faint parts of the RGB; to get enough reference stars, 
however, we had to include the upper part of the MS. (3) We 
draw a fiducial line along the MS, SGB, and RGB using refer- 
ence stars. (4) Our actual DR correction makes use of the fact 
that DR varies systematically with position in the field. What 
we do is to correct the color of each star by the median redden- 
ing of the nearest «, reference stars (see below). (5) We rotate 
the corrected values (i.e., corrected x and unchanged y) back 
into the CMD frame. 

A final step in each iteration was to redraw the fiducial se- 
quence in reaction to the changes in the choice of reference 
stars. (In practice, changes after the first iteration are trivial, 
though.) 

For any CMD that we need to correct, we choose by trial 
and error, as a compromise between a large n, for robustness 
and a small one for spatial resolution. Typical «» values are in 
the range 50-70. 

The reason we need to iterate the process is to optimize the 
choice of the set of reference stars. As this choice improves, 
the DR correction gets better, and the sequences in the cor- 
rected CMD get sharper, allowing a more perceptive choice of 
reference stars — and so on, until the choice of reference stars 
no longer improves. It is important to keep in mind, however, 
that the DR correction is always made to the original photom- 



etry, but with the new reference stars; it is the optimization of 
the reference stars that makes the method work. 

We illustrate this procedure with our DR correction of 
NGC 6441, several steps of which are shown in Figure |2] The 
first panel shows the mpguw vs. mpsgow - wifsuw CMD of the 
cluster around the SGB region, with our initial choice of ref- 
erence stars marked in red, and our initially-drawn fiducial se- 
quence in green. 

The next panel shows the result of the first application of 
the procedure; the DR correction was good enough to produce 
a significant narrowing of all the sequences in the CMD; more- 
over, the RGB is now bimodal, and the MS is showing sugges- 
tions of bimodality. This improvement allows us to choose a 
new set of reference stars (marked in red), from which we have 
removed stars of the blue MS, and a number of binaries and 
field stars. The green line is the redrawn fiducial sequence. 

As indicated, we iterated this procedure. Three iterations 
were enough to reach our final choice of reference stars and 
the consequent DR correction. (A fourth iteration yielded neg- 
ligible improvement.) The final DR-corrected CMD is shown 
in the right panel of Fig. |2] 

Our procedure also allows us to construct reddening maps, 
by averaging the DR value used to correct the stars in each 
region. (Note that we will return to this subject in Sect. 13.31 to 
present a specific validation of the NGC 6441 reddening map.) 
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Finally, we note that from this point on, all the CMDs that 
we show are corrected for differential reddening. 

2.6. Proper motions 

As noted in an earlier subsection, field-star contamination 
presents a serious problem in our analysis of the fine struc- 
ture of the CMDs. There is no problem for the ACS/HRC 
and WFPC2 fields, which are at the cluster centers, where field 
stars represent a negligible fraction. The issue arises when we 
analyze ACS/WFC data (in particulai- GO-10775) and the new 
WFC3 exposures (GO- 1 1739), which were taken off-center to 
avoid crowding. Since the GO-10775 and GO-11739 expo- 
sures overlap for both clusters, we were able to measure proper 
motions, with a baseline of at least 4 years. In addition, some 
of GO-11739 exposures of NGC 6441 were purposely taken 
a year apart, in order to measure proper motions in the par- 
allel ACS/WFC field, where it was clear that field-star con- 
tamination was very serious. This provides us with up to 5 
years of time baseline between GO-10775 and GO-11739 for 
NGC 6441 (and of course the 1 -year baseline within GO-11739 
— needed for CMDs that use only that program). 

Proper motions were measured using ACS/WFC F606W 
and F814W exposures as first epoch and WFC3/UVIS F390W 
exposures as second epoclH We did not use WFC3 IR expo- 
sures for proper motions, because of the more than twice-as- 
large pixel size of the IR detector and its severe undersampling. 

We used the catalogs created by Anderson et al. (2008il to 
define a master coordinate frame. As we did in (e.g.) Bellini et 
al. (2009a), we first defined for each cluster a reference set of 
probable members, on the basis of their position in the mpsMw 
vs. OTF606W - 'WF814W CMD. These included only bright, unsat- 
urated stars, which have smaller position errors. 

For each star we took each and every pairing of a position 
measurement from the first epoch and another from the second 
epoch, and transformed each of the two catalog positions into 
the master frame, by means of a six-parameter linear transfor- 
mation using only a local set of 50 reference stars; we then 
differenced the two epochs to derive a displacement for that 
pairing. We took a 3o--clipped median of all these displace- 
ments of the star, and assigned it a sigma that was the 68.27 
percentile of the residuals around the median. Since the refer- 
ence stars are nearly all cluster members, the displacements of 
stars are all relative to the mean motion of the cluster. 

After this first estimate of the proper motions, we excluded 
from our selection of reference stars those whose motion dif- 
fered from the mean by more than 3 sigmas, and re-derived the 
motions. After three iterations the process converged. 

Even though the mean motions of both NGC 6388 and 
NGC 6441 are similar to the motions of the field stars, the dis- 
persion of internal motion in each cluster is about five times 
smaller than that of field stars (~18 km/s vs. 100 km/s). Thus 
by selecting stars whose motion is close to the mean motion of 
the cluster, we exclude the majority of field stars. 



^The only exception is tlie proper-motion cleaning of the GO-11739 CMD of 
NGC 6441, using only F390W images, with a 1-year time baseline, 
^http : //www . astro . u£l . edu/~ata/public Jistgc/databases . html 



For NGC 6388, the overlap region can be seen in the left 
panel of Figure 1 , involving the bifurcated red and the bifur- 
cated blue footprints at the lower left. In Figure |3] we illus- 
trate our proper-motion selection. Successive panels show the 
CMD of all the stars, the separation of cluster and field based 
on the sizes of motions, the CMD of the cluster stars, and the 
CMD of the field stars. From the last panel it is clear that in 
CMDs where we do not have proper-motion separation, the 
risk of field-star contamination will be greatest for the SGB 
and MS regions. (Note: although we had three filters avail- 
able, we have chosen here to display the CMD using the more 
familiar F606W-F814W pairing; F390W will come into play 
in our ensuing discussion of the details of the sequences of the 
CMDs.) 

Figure|4]is the analog of Fig. [3] for NGC 6441, in the over- 
lap region of the same two programs (the same colors and 
shapes of footprint, near the top of the right-hand panel of 
Fig. 1). In this case there are many more stars, both because 
of a larger overlap and also because these exposures are near 
the cluster center. Also, field-star contamination is much more 
severe than for NGC 6388, making cluster-field separation for 
NGC 6441 even more important. Nevertheless, our proper- 
motion separation yields an almost clean sample of cluster 
members in panel (c), even though panel (d) shows that we 
have lost a few cluster members. 

Proper-motion-selected CMDs will be used extensively 
throughout the paper Note also that when proper-motion mea- 
surements were available we used them to select bona-fide 
cluster members as reference stars for improved DR correc- 
tions. 

3. Evidence of multiple populations 

In this section we marshal the resources of our cameras and 
filters to examine each of the clusters in detail, in order to 
uncover all the evidence of multiple populations that we can 
find. We begin with an overview of a carefully chosen array of 
CMDs of each cluster; then we take up the regions of the CMD 
in detail, one by one. 

3.1. Detailed color-magnitude diagrams of NGC 6388 

Figure |5] begins with two overall CMDs of NGC 6388. In 
panel (a) we show the most familiar form of CMD, using 
F606W and F814W from the ACS/WFC images of GO-10775; 
it extends from the tip of the RGB down to five magnitudes be- 
low the MS turn-off. (These include the stars that we showed 
in the nearly identical CMD of panel 4(a), but now we include 
the fainter stars that did not figure in the proper- motion discus- 
sion.) Panel (b) shows a CMD from F390W and F160W of the 
WFC3/UVIS and IR respectively, from GO-11739; there are 
fewer stars here, because the smaller WFC3 images are off- 
center, but this is our widest color baseline, and this CMD is 
new. 

In panels (c)-(g) we show in detail the five regions that are 
outlined in panel (a); for each detail panel, however, we use 
the combination of filters that best enhances the features of that 
part of the CMD. We will return to these in later subsections. 
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Fig. 5. — Overview of the main features in the CMDs of NGC 6388, from different combinations of magnitude and color Panel (a) 
shows a CMD from GO-10775 (the most complete CMD), while in (b) is the CMD from GO-11739 WFC3 (the CMD with the widest 
color-baseline). In (a) we have highlighted five regions, labeled (c)-(g), around main cluster features. These regions are then shown in 
panels (c) to (g) using the filter pairs that best enhance each CMD region, e.g., a double RGB in panel (c), or a double SGB in panel 
(f). 



but even now we readily see the splits of the RGB in (c) and the 
SGB in (f) as well as the complex structure of the red clump 
of the HB in (d). (It is not obvious, but we shall see that panel 
(g) shows a greater MS width than photometric errors can ac- 
count for. In any case, the corresponding panel for NGC 6441 
will show a clear split.) Panel (e) shows the extended blue 
HB, which has the anomalous blue tail of stars (highlighted by 



larger dots). In each panel the reddening vector is shown as a 
red arrow. 

In these panels we have plotted only those stars less affected 
by the severe crowding in these central fields. We omitted stars 
closer to the center than a magnitude-dependent limiting ra- 
dius, which we chose by hand in a way that excludes larger 
regions for the fainter stars. 
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NGC 6441 overview 
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Fig. 6. — Overview of the main features in the CMDs of NGC 6441, as in the previous figure for NGC 6388, except that the examples 
of enhancement are now the double RGB in panel (c) and the double MS in panel (g). 



Proper-motion selection is not used in any of the panels of 
this figure. The purpose of panels (a) and (b) is simply to show 
the whole CMD and to define panels (c)-(g). Panels (c), (d), 
and (e) show CMD regions where field-star contamination is 
minimal, while in (f) and (g) any gain in sharpness would be 
greatly outweighed by the loss in star numbers from cutting 
back to the small overlap region between the GO- 10775 and 
GO-11739 images. 



3.2. Detailed color-magnitude diagrams of NGC 6441 

Figure|6]is an overview of the CMDs of NGC 6441, similar 
toFig.|5]forNGC 6388. 

The CMD from GO- 10775 data alone (panel (a)) shows 
strong contamination from stars of the Galactic bulge, as al- 
ready mentioned. The 1-year time span of the GO- 1 1739 data 
set allowed us to remove most of the field stars from the largest- 
color-baseline CMD shown in panel (b). (This is why there is 
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less field contamination here than in Fig. [Sj^ ) Panels (a) and 
(b) here have the same scale as those of Fig. |5] so as to allow 
a direct visual comparison of the complete CMDs of these two 
clusters. (Similarly, in the following subsections we will show 
a direct comparison of the different evolutionary sequences of 
these clusters using CMDs with the same filters and axis scale.) 

Again, in the rripsuw vs. »Jf606w-'Wf814w plane we highUght 
five regions, labeled (c) to (g), where the main cluster CMD 
features are present. The CMDs zoomed around these regions 
are shown in the other panels of the figure, using the combi- 
nations of filters that best enhance these features, and taking 
advantage of proper motions to remove field stars, whenever 
useful (panels (c), (d), (f), and (g)). 

We did not use proper-motion selections for panel (e). First, 
the blue HB is only marginally contaminated by field stars. 
Second, because proper motions are confined to the over- 
lap area between GO-10775 and GO-11739, we would have 
ended up with only half as many stars, in a region of the 
CMD that is thinly populated. We need also to note that 
for panel (d), for the RC of the HB, we had to use here the 
wipSMw vs. OTF606W WJF8I4W CMD instead of the mF435w vs. 
'WF435W - 'npsssw CMD that we used for NGC 6388, because 
for NGC 6441 the equivalent WFPC2 F439W images do not 
go deep enough. 

Unlike NGC 6388, the RGB of NGC 6441 clearly spUts at 
its base when using the F390W filter (panel (c)). Furthermore, 
for the first time we can see that the MS of NGC 6441 is split 
into two distinct sequences (panel (g), more details later). 

In the following subsections we will compare and analyze 
each evolutionary stage of the two clusters. We want to em- 
phasize once again that all the CMDs presented in this work 
are corrected for differential reddening effects. 

3.3. The Main Sequences 

Figure |7] shows the CMDs of NGC 6388 (on the left) and 
NGC 6441 (on the right) around the MS region in the mpagow 
vs. mF39ow - 'MpSMw plane, adopting the same axis scale. The 
reddening vectors are also displayed. It is clear from the figure 
that the MS morphology of these "twin" clusters is quite dif- 
ferent. The MS of NGC 6388 is wider than expected from ob- 
servational errors alone (more in the following), but the CMDs 
made using all combinations of filters at our disposal fail to 
provide any hints of a MS split. On the other hand, there seems 
to be no doubt that the MS of NGC 6441 is bimodal. 

Because of the relatively high reddening values of the two 
clusters, one might suspect that the spread/split in the MS of 
the two clusters is the result of badly corrected differential red- 
dening. We note, however, that the reddening vector runs par- 
allel to the MS, so that reddening cannot affect the width of 
the MS, nor should DR correction affect it either This near- 
parallelism of the reddening vector to the MS should be kept 
in mind, as it applies to nearly all CMDs. It will be seen, 
moreover, that a similar parallelism holds for the reddening 
vector and two-color diagrams of MS stars for all color combi- 
nations. Indeed, if we plot an uncorrected CMD of a region of 
NGC 6441 small enough that differential reddening cannot be 



an issue, the MS split is visible. 

To further verify that we are not being deceived by seri- 
ous systematics in our DR corrected photometry, we performed 
two specific tests. The left panel of Fig.[8]shows the mpsgow vs. 
WJF390W - WJF814W CMD of NGC 6441 around its MS and TO 
regions, for all stars in common between GO- 1 1739 and GO- 
10775 FoVs. (We chose these two filters because they are the 
only ones that clearly split the MS; we are confined to the stars 
in the overlap region because the two filters come from those 
two different programs.) 

The inset shows the color-coded DR map for all the stars, 
in the coordinate reference system of the GO-10775 data set, 
with map coordinates in ACS/WFC pixels. The cluster cen- 
ter is marked by a "+". Each point corresponds to a star, and 
is colored according to the DR correction found for it. The 
maximum extent of the DR corrections is 0.86 mag along the 
direction of the reddening vector 

Since the DR correction map shows a vertical gradient, for 
the first test (TEST A) we split the FoV into upper and lower 
subsets containing the same numbers of stars, and performed 
the DR correction procedure on each of the two independent 
subsets, as a test of whether the apparent breadth of the MS 
couldbe just high/low reddening. The resulting CMDs, in pan- 
els (b) and (c), each show the MS split of panel (a) — some- 
what more clearly in (c), because the stars of panel (b) suffer 
photometrically from the crowding in the cluster center In the 
less perturbed photometry of panel (c) one can even discern a 
dominant red component and a bluer tail in the colors. But the 
main point is that Test A totally refutes the idea that the bluer 
and redder parts of the MS spread could simply be stars with 
greater and less reddening. 

To investigate whether the MS spread/split might in fact re- 
sult directly from crowding in the cluster center, we performed 
TEST B, in which the two halves of the sample were an inner 
one and an outer one, as shown in panels (c) and (d). The MS 
width is still evident in both samples, and the outer sample, in 
panel (d), if anything shows the split even more strongly than 
panel (c), above it. 

These tests show that the MS spread of NGC 6441 is real. 
Moreover, when only less-crowded, well-measured stars are 
selected (e.g., right panel of Fig. |7] but see also panel (a) of 
Fig.[8]l, the MS splits into two branches, the redder one being 
the more populated. 

No MS split was visible in the case of NGC 6388, but 
the tests just performed can be taken as evidence that the MS 
broadening observed there is also real, rather than an artifact 
of differential reddening. 

3.3.1. The MS of NGC 6388 

As shown in the study by Anderson et al. (2009) of the 
MS broadening in 47 Tuc, and by Milone et al. (2010b) 
for NGC 6752, an effective way of testing for a true CMD- 
sequence broadening is to divide the images of the field of stars 
into two independent halves. If the broadening is intrinsic, the 
stars that have redder (or bluer) colors in a CMD made from 
one haff of the images will have redder (or bluer) colors in the 
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Fig. 7. — mF39ow - '«F8i4w CMD of the upper-MS and SGB regions of NGC 6388 (left) and NGC 6441 (right), on identical scales, 
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CMD from the other half. But if the broadening is due only to 
measuring errors, a star that is redder in the first half of the im- 



ages will have an equal chance of being redder or bluer in the 
second half, and the bluer stars of the first subset are equally 
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the text, in which we divide the images into two sets and show that the stars show similar color residuals in the two independent data 
sets. 




1 1 1 M 1 1 m; 1 1-1 M 1 1 1 M 1 

(a) i],. 


1 M M 1 1 1 M l r.M 1 1 1 M 1 1 

-(b) 

r 

* ■" . 

_:U "■. 

.7^ '■ 
. il* ■ 


4" " 1 " " 1 " "1 ' " T 

;(c) .■ ; 
\ / ] 

i- 11 11 1 11 11 1 11 11 1 11 11 f 


1 1 M 1 1 1 1 M, M M 1 M M 1 

-I'll OUT 

S .. . 

0.203" 

f ' ■ 

■:'.o.245" 
. . XV ? ■ 


■ ^0 012" 

\ 
i 

i 

---I--- 

? (7 = 

j 0.013" 

f 
1 

i 

i 0.026 

■s. 
■'i. 


4,, M 1 M M 1 , , M 1 , M 1^ 
t 11 11 1 11 11 1 11 11 1 11 11 f- 


4" " 1 " " 1 " " 1 " "t 

L X J 


■1 , , , , 1 , , , , 1 ,■, .,J , , , , 1- 


-1 . . . . 1 . r. 1 1- 


4"" 1 " " 1 " " 1 " "4 

rx: 


■■" * .0= 

. ■ '. •.■'■0.352" 
■•■ ■ 

1 1 M 1 1 1 1 M l-l-l .iJ LI M b 


0,044" 



F390W F160W 



™F390"W -^FieOW 



-1 -0.5 0.5 1 -1 -0.5 0.5 1 -1 -0 5 0.5 1 -1 -0 5 0.5 1 -1 -0.5 0.5 1 

A; A3 A, VS. A3 (Aj + A.)/2 (A,-A3)/2 



Fig. 10. — Same as Fig.|9] but using mpieow vs. mpsgow - '«fi60w. 



likely to be red or blue in the second half. Again, we note that 
the reddening vector is almost parallel to the MS (see, e.g., 
panel (a) of Fig. |5]l, and therefore DR cannot change the lo- 
cations of the stars with respect to the MS fiducial line. In the 
following, we will apply this technique, basically following the 
procedure described in Anderson et al. (2009). 

The leftmost panel of Fig.|9]shows the mpeoew vs. mpeoew - 
wipSMw CMD of NGC 6388. The second panel shows a zoom- 
in of the upper part of the MS (colored in red in the first panel) . 
The continuous curve in red marks the adopted MS ridge line. 
"We subtracted from the color of each star that of the MS ridge 
line at the same magnitude level. The rectified MS obtained 



in this way, for the two independent set of images, is shown 
in panels (a) and (b). In panel (a) the stars were divided into 
three equally-populated groups, color-coded in red, green and 
blue. In panel (b) we have kept for each star the same color 
that it had in panel (a). (We decided to divide the MS into 
three groups rather than two, as we had done in our previous 
papers, to better understand the color behavior of the stars in 
the two subset of images.) 

It is clear that the color of stars is maintained very well (with 
only a small scatter due to photometric errors) from panel (a) to 
panel (b). To emphasize this fact, we show in panel (c) the cor- 
relation between the color of each star in the two halves of the 
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Fig. 1 1 . — Multicolor analysis for the MS of NGC 6388. Panel (a): the mpieow vs. mpsgow - '«fi60w CMD. The inset shows the rectified 
CMD, with blue, black, and red stars plotted in those colors. Panels (b) through (f) show CMDs using different filter combinations; 
the stars keep the same color coding; panel (g) is the mpscjow - 'Wf606w vs. mpguw - 'wpieow two-color diagram, showing only the blue- 
and red-coded stars. The two extreme groups remain distinct in all panels. Each panel includes a reddening vector. 



data, divided in four different magnitude bins (since measuring 
errors increase at fainter magnitudes). A star redder than the 
MS ridge line in the first half of the data is also measured as 
redder in the great majority of cases, in the second half of the 
data, and vice-versa. This is the mark of a true spread in color 

In panel (d) we show the color distribution of the rectified 
MS as measured from the whole set of images, (Ai + A2)/2. 
The distribution of half the difference between the colors mea- 
sured separately from each half of the images, (Ai - A2)/2, is 
presented in panel (e). This distribution is a statistical estima- 
tor of the error in the color distribution in panel (d). 

At each magnitude level of panels (d) and (e), an estimate 
of the cr of the spread of the distribution is given (calculated 
here as the 68.27 percentile of the absolute deviation from the 
median). As one can see, the intrinsic width of the MS is about 
three times as large as what would be expected from photo- 
metric errors alone. (As a further test, artificial stars within 
the magnitude interval shown in Fig. 9, and chosen to have 
'Mpeoew - wipSMw colors along the MS of NGC 6388, gave MS 
widths ranging from 0.005 to 0.027 — only a small fraction of 
the width observed in Fig. 9.) 

A similar procedure applied to the nzF390w - nzF8i4w color 
yields a Icr width of the MS of 0.09 mag, at 1 . 1 mag below the 
turnoff (after subtracting in quadrature the photometric errors 
of ~ 0.03 mag). 

The presence of an intrinsic breadth in the MS in NGC 6388 
is further confirmed when we apply the technique of Anderson 
et al. (2009) to the GO-11739 data set (Fig.[TSll. In this case, 
the color baseline is six times larger, and the MS appears to be 
about ten times wider than what is expected from photometric 



errors alone. 

The GO-10775 and GO-11739 fields overlap marginally 
(only about a thousand stars with a measure in all the 4 avail- 
able filters), but this is enough to show, in Fig.[TTl a multicolor 
analysis, similar to what was done for d Centauri by Bellini 
et al. (2010), 47 Tuc by Milone et al. (2012b), and NGC 6397 
by Milone et al. (2012c). Moreover, here we can use proper 
motions to eliminate field stars — much more important here 
than in the sparse fields of Milone et al. 

We rectified the MS in the wzpieow vs. »zp39ow - '«fi60w 
plane, in the magnitude interval 18.75 < mpieow < 21 (see 
panel (a) and the inset of Fig. fTTT i. Note that at wpigow > 
20.5 the MS is already as wide as about 1 magnitude in the 
proper-motion-selected CMD. We defined three groups of stars 
(marked with blue, black, and red colors), with the property 
of being bluer, of average color, and redder, respectively, than 
an arbitrary line drawn through the middle color of the MS 
(not shown in the Figure). We then kept the same colors for 
the stars of these two groups in plotting different CMDs (pan- 
els (b) through (f)), made by using all other available color 
combinations. A final panel (g) shows the two-color diagram 
'MFsgow — 'WF606W VS. mp8i4w - OTpieow for the two extreme 
groups of stars. It is interesting to note that blue and red stars 
always remain on the blue and on the red side of the MS, re- 
spectively. The reddening vectors are drawn in order to empha- 
size that they are approximately, or even closely, parallel to the 
faint MS in the CMDs and, most importantly, in the two-color 
diagram. Thus differential reddening has little or no effect on 
the MS width or on the separation of the red and blue points. 

We have demonstrated that the MS of NGC 6388 shows a 
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significant color spread. This can be seen both by using opti- 
cal colors and, more importantly, by using the wide color base- 
line offered by the new GO-11739 F390W and F160W images. 
Even so, our careful investigation was unable to reveal any hint 
of MS splitting with the data sets at our disposal. 

3.3.2. The MS of NGC 6441 

In NGC 6388 we were able only to demonstrate that the 
MS is broad, without managing to reveal possible MS sub- 
populations. In NGC 6441, by contrast, the MS components 
are visibly split, and we can therefore study their behavior in 
a much more straightforward way. What is more, GO- 10775 
and GO-11739 overlap by about half the ACSAVFC field — 
much more than in the case of NGC 6388 — and this makes 
for much better statistics. What we do for NGC 6441, then, is 
to perform a multicolor analysis in order to model the differ- 
ence in the chemical content of the two MS populations, as we 
did in 47 Tuc (Milone et al. 2012b) and NGC 6397 (Milone et 
al. 2012c). 

We drew a fiducial line by hand passing in between the two 
MSs, and defined as rMS stars those within a 0.2-mag interval 
of color redward from this fiducial, and as bMS stars those 
within 0.2 mag blueward from it, in the magnitude interval 
23.5 < mF39ow < 26.5 (see panel (a)) of Fig.[T2b. 

Panels (b) show the mpguw vs. mx - mpsMw CMDs of the 
two MSs, where X is one of the other three filters (F390W, 
F606W, and F160W). In panels (c) we plotted the location of 
the ridge lines of the two MSs in the same planes as panels (b). 
We measured the color difference Acolor between the blue and 
the red ridge lines at a fixed magnitude level (mpguw — 21.5, 
black horizontal lines). Panel (d) collects these color differ- 
ences as a function of the central wavelength A of the nix band 
(open circles). 

The color difference between the two sequences depends 
strongly on the adopted color baseline. In particular, it is high- 
est in the mpjgow - 'Wf814w color, where it amounts to ~ 0.166 
mag (1.1 mag below the turnoff). This compares to 0.09 mag 
as the width of the MS of NGC 6388 and we conclude that the 
MSs of the two clusters are intrinsically different, with the MS 
of NGC 6388 being narrower than the separation of the two 
MSs of NGC 6441. We will demonstrate in Section|5]that: 

• The color spread of the MSs of the two clusters is similar 
in »iF606w - '«F8i4w, implying that the difference in He 
content among the different stellar populations in the two 
clusters must be similar (as strongly suggested by the 
morphology of their HBs); 

• Other chemical differences (which would be C, N, and 
O) must be invoked to explain the sharp difference of the 
two MSs in the mpguw vs. mpsgow — 'WpgHw CMD. 

3.4. The Sub-Giant Branches 

The top panels of Fig. [13] show the CMDs of NGC 6388 
and NGC 6441 around the SGB region in the »zf814w vs. 
wpeoew - WZF814W plane, both with the same axis scale. All 
the stars plotted come only from GO- 10775. While the SGB 



split in NGC 6388 is clear, field-star contamination prevents us 
from analyzing the SGB morphology of NGC 6441 in detail. 
Nevertheless, it is clear that the overall width of the SGB of 
NGC 6388 is decidedly larger than that of NGC 6441 . 

Thanks to the overlap between the fields of view of GO- 
10775 and GO-11739, we can take advantage of the ~5-year 
time baseline between the two and use proper motions to iso- 
late probable members of NGC 6441; panel (c) shows the re- 
sulting CMD. Now that field-star contamination is no longer 
an issue, there seems to be a second, fainter and less populated 
SGB even in NGC 6441 (as could already be suspected in panel 
(f) of Fig.©. 

Panel (d) is a replica of panel (c) on which we have drawn 
(in red) the ridge line of the MS, SGB, and RGB as follows. 
We divided the cluster stars in the figure into 0.15-mag bins 
(0.4 mag for RGB stars) and derived 3cr-clipped median val- 
ues of magnitude and color for each bin. We selected only the 
brighter SGB for this purpose, while stars that seemed to be- 
long to a fainter SGB have been colored in red. Using the same 
technique, we derived in a similar way for NGC 6388 the ridge 
lines for the MS, SGB (both components), and RGB, and drew 
them in blue. We manually shifted the main-component ridge 
line of NGC 6388 to best overlap the MS and turn-off regions 
of NGC 6441. 

The SGB split of NGC 6441 is much less than that of 
NGC 6388. This means that, if a difference in C+N+O is 
the cause for the SGB split in NGC 6441 and NGC 6388 (as 
suggested for the case of NGC 1851 by, e.g., Cassisi et al. 
2008), that difference must be smaller for NGC 6441 than for 
NGC 6388. 

We performed a multicolor analysis on the SGB region of 
the clusters. A close-up view of the SGB of NGC 6388 is 
shown in panel (a) of Fig. [141 for proper-motion-selected stars 
in common between GO-10775 and GO-1 1739. The blue and 
red colors are used to define the bright and the faint SGB com- 
ponents, hereafter bSGB and fSGB, respectively (i.e., in the re- 
maining panels each star will keep the color that it is assigned 
here). Panels (b) through (f) show the other five CMDs made 
with the remaining color combinations at our disposal (six to- 
tal combinations with four filters). The reddening vector is 
displayed in all panels. Bright and faint SGB stars are color- 
coded as in panel (a). In all panels the bSGB is brighter than 
the fSGB, and the two components remain well separated. 

Similarly, Fig.[T5]shows the multicolor analysis for the SGB 
of NGC 6441. The bright and faint SGB components for 
NGC 6441 are selected in the mpeoew vs. mF606w-'WF8i4w CMD 
(panel (a)), and again we plot only proper-motion-selected 
stars. 

Unlike the case of NGC 6388, the two SGB components of 
NGC 6441 do not remain distinct in all CMDs, but the bSGB 
stars again stay, on average, brighter than the fSGB stars. The 
separation of the two groups, however, is much less obvious 
than for NGC 6388. 

A deeper analysis of the SGBs of NGC 6388 and other glob- 
ular clusters is presented in a separate paper (Piotto et al. 2012). 
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Fig. 12. — Multicolor analysis of the MS of NGC 6441. The two MS populations are selected in panel (a). Panels (b) and (c) show 
the behavior of the sequences and their fiducials in the planes mpgww vs. nix - 'Wf814w, where X is one of the other filters. The color 
differences between the two sequences at fixed mpsMw = 21.5 are plotted in panel (d). 



3.5. The Red-Giant Branches 

Both NGC 6388 and NGC 6441 have a split RGB, but the 
splits are different, as shown in Figure [16] in which we see 
CMDs of the two clusters side by side, with two different color 
baselines, top and bottom. Unfortunately, in CMDs drawn 
from quite disparate sources the numbers of stars are not com- 
parable; moreover, the presence of the asymptotic giant branch 
(AGB) next to the RGB confuses the picture further Neverthe- 
less, we shall show how to use these CMDs to trace where the 
brighter and fainter RGB sequences (bRGB and fRGB) run in 



each cluster 

In panel (c) of Fig.[T6]a spHt in the RGB of NGC 6388 is ev- 
ident for stars brighter than mpsssw = 16.2. We can exclude the 
possibility that the brighter RGB sequence might be composed 
of AGB stars, since the latter are bluer and well separated from 
the RGB. Panel (a) seems to be suggesting a separation for the 
entire length of the RGB, but it is hard to draw definite conclu- 
sions from so few stars. 

For NGC 6441 panel (b) shows a clear split in the faint part 
of the RGB, but the split becomes less convincing at brighter 
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Fig. 13.— (a) and (b): The SGB region of the CMDs of NGC 6388 (left) and NGC 6441 (right), showing all stars from GO-10775. 
(c): Proper-motion-selected members of NGC 6441, from the overlap region of GO-10775 and GO-11739. Panel (d) is a replica of 
panel (c), but with the ridge lines of the MS, the bright part of the SGB, and the lower RGB of NGC 6441 marked in red, and those 
of NGC 6388 over-plotted in blue — the latter shifted to fit the stars of NGC 6441. The stars that are colored red are the fSGB that is 
discussed in the text. 



magnitudes. Panel (d) offers weak confirmation of a split at the 
bright end, but the star numbers are just too small to be at all 
sure. 

Interestingly, RGB splits are usually found only in ultravio- 
let photometry (or else in Stromgren filters that are sensitive to 
molecular bands see, e.g., Gratton et al. 2012, and references 
therein, for a review); panels (a) and (b) are in agreement with 
this. By contrast, the cluster M 4, for example, has a very nar- 



row RGB in the B - / color (Marino et al. 2008), but is broad, 
or even split, in U - B. Thus the apparent split of the RGB of 
NGC 6388 in an optical color (panel (c)) comes as a surprise, 
and may be suggesting an additional population phenomenon 
in NGC 6388 and NGC 6441, different from that at work in 
M 4. Unfortunately, all the stars for which an abundance anal- 
ysis have been performed (e.g., Carretta et al. (2007), lie at 
radial distances outside our FoV, and therefore we cannot link 
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Fig. 14. — (a): mpsuw vs. mpeoew - '«f814w CMD of the SGB region of NGC 6388 for proper-motion-selected stars in the overlap 
between GO-10775 and GO-11739, with bSGB and fSGB stars colored azure and red, respectively. The remaining panels show the 
other color combinations that are available, with each star retaining the color that it was assigned in panel (a). The bSGB and fSGB 
stars remain well separated in all panels. Reddening vectors are also shown. 



a given RGB branch to a particular Na-O family. 

In conclusion, the RGBs of these two clusters show evi- 
dence of splits, but firmer conclusions await richer samples of 
stars, or a wider choice of filters, or both. 

3.6. The Horizontal Branches 

We will focus our HB analysis only on the red clump (RC) 
region, and not on the extended HB (which is discussed in Sec- 
tion |52] instead, in conjunction with the helium abundance). 

Milone et al. (2012b) showed that UV two-color diagrams 
allow a clear separation of the two main populations of 47 Tuc 
in all evolutionary sequences, and we can profitably use a 
similar technique on the RC of NGC 6388. Here we chose 
the small but high-resolution ACS/HRC exposures in F330W, 
F555W, and F8 14W. This data set was the best choice we could 
make. We tried adding F435W from WFPC2 SNAP-9821 for 
NGC 6388, but it turned out not to add any useful information. 
We could not combine ACS/HRC with ACSAVFC from GO- 
10775, because of the severe crowding in the latter Had we 
tried instead to use WFC3 F390W exposures, combined with 
GO-10775, the small spatial overlap would have given too few 
stars for adequate statistics. 

The top row of panels of Fig [IT] shows the three possible 



two-color diagrams from the HRC data set. Panels (dl), (d2), 
(el), (e2), and (fl), (f2) show the CMDs of the RC for different 
color-magnitude combinations. 

The choice of red-clump stars that we made is shown by the 
box in panel (fl). In panel (a) we drew by hand a line that 
separates the stars into a main and a secondary group (mRC 
and sRC). The groups are colored red and blue, respectively, 
and each star retains its color from this panel throughout the 
rest of the figure. 

Panels (a) and (c) clearly show a double-horned RC, a fea- 
ture that is less obvious in panel (b). We found the sRC to be 
fainter overall than the mRC (with perhaps the exception of 
panel (el)). Moreover, sRC stars are bluer than mRC stars in 
panels (dl) and (el), but redder in (fl) and (f2). 

Anomalies in the RC morphology have been reported also 
in the CMDs of NGC 6440 and NGC 6569 by Mauro et al. 
(2012) and, on a different scale, in the CMD of Terzan 5 by 
Ferraro et al. (2009), both using IR filters. Since here we in- 
stead see the RC split in a two-color diagrams employing UV 
filters, we think that the cause of the RC split in NGC 6388 
is more similar to that of, e.g., 47 Tuc (Milone et al. 2012b), 
rather than those seen in IR CMDs. 

We performed the same multi-color analysis also for the RC 
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Fig. 15. — Same as Fig. [141 but for the SGB of NGC 6441. (There are more stars because the overlap between GO-10775 and in 
GO-1 1739 is larger for NGC 6441.) 



of NGC 6441. To better compare the results between the two 
clusters, we made use of observations taken through filters as 
similar as possible to those used for NGC 6388. This choice 
implied the use of the WFPC2 observations from GO-5667 and 
GO-8718, which include data collected with the F336W filter 
instead of the ACS/HRC F330W one. Because of this, results 
on the RC analysis of NGC 6441 are not one-to-one compa- 
rable with those of NGC 6388, but they can still tell us some- 
thing about similarities or differences between the RC of the 
two clusters. 

Figure [TSlis the equivalent of Fig. [TTlfor NGC 6441. Red- 
clump stars are again selected in the gray-line box in panel (f 1 ), 
and stars of the mRC and sRC are again delineated in panel (a) 
and carry the same colors throughout. 

The first difference from the RC of NGC 6388 is that here 
none of the two-color diagrams has a double-horned shape. In 
all the CMDs, stars of the mRC and sRC are better mixed than 
in NGC 6388. In panels (dl) and (d2) sRC stars appear to be 
bluer (and fainter) than mRC stars. This mRC/sRC dichotomy 
seems to carry over, although in a less evident way, into panels 
(el) and (e2) as well. 

The HE regions of these two clusters are at first glance so 
similar that many investigations (e.g.. Rich et al. 1997) lumped 
the two clusters together; nevertheless, when we look closely 
the two HBs (in particular their red clumps) are rather dif- 
ferent. [Note, however, that the double-homed red clump in 



NGC 6388 is visible only when the ACS/HRC F330W filter is 
used. Although the WFPC2 filter F336W throughput is quite 
similar to that of the ACS/HRC F330W, small differences in fil- 
ter response may produce quite large effects on the CMD (see, 
e.g., Momany et al. 2003).] We will resume the discussion of 
the HBs in Section|5] 

4. Radial distributions 

NGC 6388 and NGC 6441 have both been revealed to host 
multiple stellar populations. In this section we will study the 
radial distributions of these subpopulations, for different CMD 
regions. 

4.1. NGC 6388 

4.1.1. The Sub-Giant Branch 

We started with the double SGB of NGC 6388. We chose 
to analyze the mpeoew vs. wpeoew - 'Mfshw CMD of GO-10775 
because of the wide separation between the two SGBs in this 
CMD, and the large number of available stars. 

Figure [19] illustrates our analysis of the radial distributions 
of the SGB components. To rectify the SGB sequence we used 
a procedure developed by Milone et al. (2009), which we de- 
scribe here in some detail, as it is rather complicated. First we 
selected by hand two points (Pif, P2f) on the faint SGB (fSGB) 
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Fig. 16. — RGBs, plotted on the same scale to facilitate comparison. NGC 6388 is on the left and NGC 6441 on the right. Top and 
bottom show different color baselines. (F435W and F439W are equivalent filters in two different cameras.) Because the CMDs come 
from overlap fields of different size, NGC 6388 has fewer stars than NGC 6441 in the top row and more in the bottom row. 



and two points (Pib,P2b) on the bright SGB (bSGB). These 
points define the two Unes in panel (a) of Fig. [191 We then 
linearly transformed the CMD into a reference frame (which 
we will refer to simply as 'abscissa', 'ordinate') in which the 
origin corresponds to Pib and the coordinates of Pif, P2b and 
Pif are (1,0), (0,1) and (1,1), respectively. (The transformation 
equation is given in the Appendix of Milone et al. 2009.) This 
transformation equalizes the separation of the sequences at the 
two ends of the interval, but they are curved in between the 
ends. (We do not show this intermediate frame.) 



We then used the following procedure to effect a straight- 
ening that was a compromise for the two sequences. First we 
fitted a fiducial line to each of the two sequences, as follows: 
We divided the bSGB into 'ordinate' intervals of 0.025 and 
the fSGB into intervals of 0.05. For each sequence, within 
each of these intervals we computed 2.5-cr-clipped median 'ab- 
scissa' and 'ordinate' values, and then put a smooth spline 
curve through those points. Then as the line to use for rec- 
tification we took a line midway between the smooth spline 
curves of the two sequences. 
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Fig. 17. — Two-color diagrams (top panels) and CMDs (middle and bottom panels) of the red clump of NGC 6388. The RC sample is 
defined in panel (fl). The selection of the main and secondary RC components (mRC and sRC, respectively in red and blue) is shown 
in panel (a). 



The rectification itself consisted of subtracting from the ' ab- 
scissa' of each star the 'abscissa' of this line at the same 'ordi- 
nate' level. Panel (b) shows the two SGBs rectified in this way. 
For simplicity, we will refer to the newly obtained A'abscissa' 
as just 'abscissa'. We also show in panel (b) the transformed 
locations of Pit, P2f, Pib, and P2b. 

Stars with 'abscissa' between -1.6 and 1.6 (gray vertical 
lines) and 'ordinate' between 0.1 and 0.85 (in red) were se- 
lected for the radial-distribution analysis. 

We excluded from the analysis the innermost 25 arcsec be- 
cause of high crowding, and divided the remainder of the GO- 
10775 field into three radial intervals (marked 1, 2 and 3 in 
panel (c)) in such a way as to have the same number of stars in 
each interval. The excluded innermost 25 arcsec are colored in 
azure. 

For each radial interval (panels (dl), (d2), and (d3)) we ex- 
tracted the 'abscissa' distribution of the stars within the limits, 
and did a simultaneous fit of a pair of Gaussians to it, as done 
in Bellini et al. (2009b). The fit was limited to the 'abscissa' 
values between -1.6 and 1.6. The fit is marked with a black 
curve, while the individual Gaussians for the bright and the 



faint SGB are colored in blue and red, respectively. In the top- 
left corner of each of these panels we reported the areas of each 
of the two Gaussians (A^bSGB and A^fsoB, for the bright and the 
faint component, respectively). 

After this somewhat complicated fitting procedure, it is not 
immediately obvious how to assign error bars to the numerical 
ratios of the two components, but the answer becomes clear 
when we step back and consider what we have done. Our basic 
aim has been to determine, at each distance from the cluster 
center, what fraction of the population belongs to each compo- 
nent; the fitting of the dual Gaussians was only a device for 
estimating how many stars belong to each population com- 
ponent. The statistics of the problem are therefore binomial: 
if a fraction pfscB of the stars actually belongs to the fainter 
component, random samplings from this population will have 
a binomial distribution, with 0-^,5^5 — ^JNp(l - p), where 
N = MsGB + M)SGB is the sample size. For each observed case, 
we use as the estimator of the fraction p the quantity pfscB = 
N(scb/N. Then crp^^^^ = o-Ntsaa/^ = V/^fSGB(l - PfSGB)/N; 
this is our best estimate of the error bar of p. 

That is what we have done. One should bear in mind, how- 
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ever, that in so far as our estimation methods might be some- 
what imperfect, the true uncertainty could be somewhat larger 
than this error bar (We prefer to make this simple statement, 
rather than attempt to doctor the statistics.) 

In panel (e) we summarize the results of the analysis: values 
of PfSGB are plotted as a function of the median radial distance 
of the stars in each radial interval, with errors. The horizontal 
lines show the radial range to which each point applies. Verti- 
cal dashed lines mark the core radius and the half-light ra- 
dius rji. We ran artificial-star experiments to check for possible 
influence of incompleteness on the radial gradients. We found 
that in the radial interval covered by our images the complete- 
ness is always higher than 84%, with marginal diff'erences (at 
most 2%) between bSGB and fSGB stars. Therefore incom- 
pleteness cannot significantly change the results displayed in 
Fig. 19. 

We notice some marginal evidence of a gradient, with fSGB 
stars more concentrated than bSGB ones. 

4.1.2. The Red-Giant Branch 

To study the radial distributions of flie RGB of NGC 6388 
we chose the ACS/WFC F435W-F555W CMD of SNAP- 



9821, the data set that was able to split the RGB in Fig. [16] 

In panel (a) of Figure|20]we see the bRGB and fRGB (as we 
again call the two components) well separated at the brighter 
end but badly blended at fainter magnitudes — a situation for 
which our dual-Gaussian technique is well suited. 

We began with a transformation and rectification procedure 
like the one that we described for the SGB. This case was a 
little easier, because the two population branches turned out to 
be parallel after the transformation, so that fewer steps were 
required for the rectification. 

A glance at panel (b) reveals a serious problem, however: 
the encroachment of the AGB on the bRGB. (We purposely 
show a wide range of 'abscissa' so as to expose this issue.) 
Because of the AGB, the bluest range of 'abscissa' will have 
to be "out of bounds" for our fitting of the Gaussians, with a 
cut-off' not far to the left of the peak of the Gaussian that fits 
the bSGB. The placement of the cutoff is then a delicate ques- 
tion, to answer which we invent a special procedure. Without 
any question of radial binning yet, we analyze all the stars that 
are between 'ordinate' values and 1 in panel (b). We make 
use of the fact that the number of AGB stars is a very differ- 
ent function of 'ordinate' from the number of RGB stars. As 
a criterion for absence of AGB interference, we divide the 'or- 
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Fig. 19. — (a) Blow-up of a region of the SGB of NGC 6388, with the four points and two lines marked that are used in the trans- 
formation to the coordinate system of the next panel. (See text for details.) (b) Rectification of the SGB region that was shown in 
(a). The horizontal black lines delimit the 'ordinate' range used in the analysis, and the vertical gray lines the 'abscissa' range, (c) 
The parallelogram is our field. The azure circle is the excluded region in the cluster center; the pink area is the working region that 
we divided into the three subregions that are marked, (d) Dual-Gaussian fits, in the 'abscissa' interval from -1.6 to 1.6, of the color 
distribution of stars in each radial subregion; 'abscissa' is the same scale as in (b). Vertical gray lines show the range that was fitted. 
Black lines are the fits, red and blue lines the individual Gaussians. (e) Relative proportions of the two populations, expressed as the 
fraction of the total that is contributed by fSGB. Horizontal lines show radial extent of each subregion; error bars come from the sample 
size and the binomial distribution. Dotted Unes mark the core radius and half-fight radius. 



dinate' range from to 1 into three separate parts, and require 
that the ratio A^fRCB/MiRGB, as derived from our dual-Gaussian 
fitting with the lowest values of 'abscissa' excluded, come out 
the same for all three ranges of 'ordinate'. We carry out this 
test (which is not illustrated here) for various placements of 
the cut-off, and choose the cut-off that is farthest to the left but 
stUl satisfies the criterion. That value is the vertical gray line 
in panel (b). 

(A thoughtful reader may ask how we can be sure that 
the two populations represented by the bRGB and the fRGB 
should have the same distribution of stars with 'ordinate' value. 
The answer is that this distribution is set by gross rules of stel- 
lar structure, as a star moves up the RGB; these rules are quite 
insensitive to the small differences that separate these two pop- 



ulations. AGB evolution is governed by quite different physics, 
however, and leads to a very different distribution of 'ordinate' 
values.) 

With this problem disposed of, we can proceed to the selec- 
tion of radial bins. Crowding is not an issue for these bright 
stars; we can study them even in the very core of the cluster. 
Farther from the center, however, the star numbers become too 
small. In panel (c) the azure outline is our SNAP-9821 field, 
and the region that we studied is in pink. 

We had found the best 'abcissa' interval to have its left limit 
at -0.6, which is marked with a gray vertical line in panels 
(dl), (d2), and (d3). (We show the full 'abscissa' range of 
panel (b), however) The best-fit dual Gaussians are shown in 
each of the panels (d), in blue and red for the bright and the 
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Fig. 20. — (a) Blow-up of a region of the RGB of NGC 6388, with the four points and two lines marked that are used in the transforma- 
tion to the coordinate system of panel (b). The analysis makes use only of stars that are to the right of the gray vertical line and between 
the two red horizontal lines in panel (b). (c) The irregular outline is the total field of view of SNAP-9821, which we used in its entirety, 
because of the paucity of stars. The three radial zones contain equal numbers of stars, (d) Dual-Gaussian fits of color distribution of 
stars in each radial subregion; 'abscissa' is the same scale as in (b). Black lines are the fits, red and blue lines the individual Gaussians. 
(e) Relative proportions of the two populations, expressed as the fraction of the total that is contributed by fRGB. Horizontal lines 
show radial extent of each subregion; error bars come from the sample size and the binomial distribution. Dotted lines mark the core 
radius and half-light radius. 



faint RGB component, respectively. We also show the areas of 
the Gaussians (i.e., star numbers) in the top-left corner 

Finally, in panel (e) we report ptHGB - A^fRCB/CMnGB + 
MjRGb) as a function of the median distance of the stars from 
the cluster center in each radial interval, with error bars c alcu- 
lated as explained in the previous subsection. 

There is no evidence for a radial gradient, though it must be 
noted that the error bars are much larger than in the previous 
case, because of the smaller number statistics. 

4.1.3. A possible relation between the radial distributions of 
the SGB and the RGB 

The radial distribution of pfsoB extends from the half-light 
radius, at Q'.5, out to ~ 2'4; our results give marginal evidence 
of a decrease (at ~2cr). Our RGB analysis is extended to larger 



radii r < 2'65, but with a much smaller number of stars and, 
for what it its worth, no evidence for a radial gradient in the 
population fractions. In both cases one component (bSGB or 
bRGB) is significantly more numerous than the other It is 
tempting to consider the bRGB as the progeny of the bSGB, 
although ^fRGB is larger than the inward extrapolation of pfscB ■ 
We will discuss possible implications of these results in the 
final section. 

4.2. NGC 6441 

4.2.1. The Main Sequence 

We have seen that NGC 6441 hosts two distinct stellar pop- 
ulations, clearly detectable both on the MS and on the RGB 
(whereas for NGC 6388 the separations were on the SGB and 
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Fig. 21. — (a) Blow-up of a region of the MS of NGC 6441. For rectification we used the average MS ridge line; the result is shown in 
panel (b). The stars selected for the analysis are in the region bounded by the red horizontal lines and the gray vertical lines, (c) The 
footprints of the GO-10775 ACS/WFC images (large parallelogram) and the GO-11739 WFC3/UVIS images (small parallelogram); 
the pink area is our working region, with our three radial zones marked, (d) Dual-Gaussian fits of color distribution of stars in each 
radial zone; Acolor has the same scale as (b). Black lines are fits, red and blue lines the individual Gaussians. (e) In black, the fraction 
of the MS total that is bMS, for the three radial intervals shown in panel (c). Horizontal lines show radial extent of each subregion; 
error bars come from the sample size and the binomial distribution. The large number of stars allowed a finer radial sampling as an 
alternative, using five subregions instead of three. Those results are shown in red. Dotted lines mark the core radius and half-light 
radius. 



RGB). 

To measure the population fraction of the two MSs we again 
used the dual-Gaussian fitting technique. Panel (a) of Fig. 1211 
shows a close-up of the MS of NGC 6441 in the mp39ow vs. 
wifsww ~'Wpgi4w CMD, chosen because this combination of fil- 
ters enhances the split of the MSs. We computed the MS ridge 
Une as done in subsection 14.1.1 1 for the SGBs of NGC 6388, 
then we subtracted from the color of each star the color of 
the ridge line at the same magnitude as the star The rectified 
MS (Acolor) is plotted in panel (b), in the magnitude interval 
23.05 < mF39ow < 27.05. Stars with 24 < mpjgow < 26 (red 
horizontal lines) and Acolor between -0. 19 and 0. 19 (gray ver- 
tical lines) were used for the radial-distribution analysis. 

Panel (c) shows the footprints of GO-10775 and GO- 1 1739 
(larger and smaller parallelograms, respectively). The overlap 
region used for the analysis is colored pink. We excluded the 



centermost 50 arcsec because of crowding problems. As done 
for the other radial-distribution analyses, we defined three ra- 
dial bins, each containing the same number of selected stars. 

Panels (dl), (d2), and (d3) show the histograms of the 
Acolor distribution within each radial interval. The best-fit dual 
Gaussians are shown in blue and red for the bMS and the rMS, 
respectively. The areas of the Gaussians are also shown, as 
done in Figs.[T9]and|20l in the top-left of each (d) panel. The 
'abscissa' interval used to fit the dual Gaussian is marked with 
gray vertical lines. 

For each radial interval we derived the values of pbMS = 
A^bMS / (M)MS + Ni-Ms), with errors, which are in black in panel 
(e). Motivated by the fact that there are more than 3000 stars in 
each radial interval, we decided also to derive a finer sampling 
of the radial distribution of /3bMS- Therefore, we divided the 
pink area in panel (c) into five equally populated radial inter- 
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Fig. 22.— As in Fig. El but for the RGB sub-populations of NGC 6441. 

the Acolor distributions for the three radial intervals, together 
with the best-fit Gaussians and their areas. In blue we mark 
the bRGB component, in red the fRGB one, and the combined 
curve is in black. 

For each radial interval, we computed pbRGB = MjRGb / (M)RGb+ 
Mrgb ) along with its error, and plotted these as a function of 
median radial distance, in panel (e). 

As for the MS, there is some evidence of a radial gradi- 
ent, with the bRGB more concentrated than the rRGB. Analo- 
gously with the SGB and RGB of NGC 6388, the similarity in 
radial gradients here is consistent with the populations of the 
RGB of NGC 6441 being the progeny of those of the MS - 
although here again there is a disagreement in the mean num- 
bers. 

4.2.3. A possible relation between the radial distributions of 
the SGB and the RGB 

In summary, our values of pbMS extend from ~0f8 (inside 
of which the crowding is too great) out to ~2'4; in this range, 
the number of bMS stars relative to the rMS shows a sig- 
nificant outward decrease. Similarly, our study of the RGB 
showed an outward decrease of the number of bRGB stars rel- 
ative to fRGB. However, as in the case of SGB and RGB stars 



vals. For each of them we fitted a dual Gaussian (not shown) 
and computed ;3bMS values, which are reported, in red with er- 
rorbars, in panel (e). There is strong evidence of a radial gra- 
dient, with the bMS more concentrated than the rMS. 

4.2.2. The Red Giant Branch 

Finally, we analyzed the radial distribution of the two RGB 
components of NGC 6441. A close-up view of the RGBs 
in the mpsgow vs. mpsgow - wJf606w CMD is shown in panel 

(a) of Fig. [22l To rectify the sequences, we derived a fidu- 
cial line for the whole RGB, and subtracted from the color of 
each star the color of the fiducial line at the same magnitude 
level. The rectified RGB is shown in panel (b). For the radial- 
distribution analysis we selected stars in the magnitude range 
22.52 < otf39()w ^ 20.84 (red horizontal lines in panel (a)), and 
between Acolor -0.13 and 0.13 (gray vertical lines in panel 

(b) ). 

Since we are using the same data set as for the MS of 
NGC 6441, panel (c) shows the same footprints as Fig.l2ll:. We 
excluded the centermost 11.1 arcsec because of crowding, and 
defined three radial intervals, each containing the same number 
of stars. 

In panels (dl), (d2), and (d3) we show the histograms of 
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in NGC 6388, here too there is an inconsistency between the 
overall levels of ^bMS and i%rcb- We will discuss the implica- 
tions of these results in the final section. 



Table 3 

Results of Radial Distribution Analyses for NGC 6388 



SGB 
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''mcd (') 


N 


PfSGB 
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A' PfRGB 
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0.711 
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0.269 


0.015 


0.195 
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0.043 
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0.014 


0.529 


144 0.361 


0.040 


1.661 


901 


0.228 


0.014 


1.141 


137 0.358 


0.041 



4.3. Population ratios and their radial gradients 

In this section we summarize the results illustrated in Fig- 
ures 19-22, concerning the number ratios between the two 
populations, and the radial gradients of the ratios. In the outer 
regions of NGC 6441 (outside ~0.'8), the blue MS represents 
~l/3 of the total population (/3bMS~ 0.33), with a small ra- 
dial gradient; closer to the center the crowding prevents us 
from distinguishing the two MSs. For the RGB we can reach 
the cluster center, where the fRGB actually predominates by a 
little, with a radial gradient bringing /?fRGB below 0.5 farther 
out. Unfortunately, radial gradients make it difficult to com- 
pare population ratios for different evolutionary sequences, but 
it does look, at first sight, as if: (i) the CMD locus of the bMS 
population crosses through the locus of the rMS population; 
(ii) the two SGBs overlap, and it is difficult to say whether the 
bMS population joins the bRGB or the fRGB. For NGC 6388 
we could not split the MS, so the population ratio and its radial 
trend were derived only for the SGB and the RGB — com- 
plicating comparisons, by leaving us with MS and RGB for 
NGC 6441, but SGB and RGB for NGC 6388. 

The results for the two clusters are summarized in Tables 
3 and 4, respectively, which give the median radius of each 
zone, the number of stars in the zone, the fraction that belong 
to the fainter component (or in the case of the MS, to the bluer 
branch), and the statistical uncertainty of the fraction. Alter- 
natively, the same results can be viewed graphically in the (e) 
panels of Figures 19-22. 

The overall tendency is to have one sequence more con- 
centrated towards the cluster center, although this is not al- 
ways statistically significant. In particular, for NGC 6441 both 
bMS and bRGB are more centrally concentrated, and the latter 
can be considered the progeny of the former For NGC 6388, 
instead, while the fSGB population is more centrally concen- 
trated, the distribution of the two populations in the RGB re- 
gion is found to be flat. Based solely on star counts, we can 
tentatively consider bRGB as the progeny of bSGB, each be- 
ing the more populated branch. 

In the study of the population ratios listed in Tables 3 and 
4, it is important to judge whether apparent radial gradients 
are real or not. To test this, we calculated a least-squares line 
through each set of points. With at most five observed points, 
any formal error for the slope of the line would be meaning- 
less, so we instead made a Monte Carlo simulation. In place 
of each observed point we took a Gaussian with mean equal to 
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Results of Radial Distribution Analyses for NGC 6441 
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Table 5 

Statistical Analysis of the Radial Gradients 



Fig. 
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0.021 
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20 
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21 
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21 
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5 
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22 
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Table 6 
multipopulation summary 



sequence 


NGC 6388 


NGC 6441 


MS 


broadened 


split 


SGB 


split 


hints of a split 


RGB 


split 


split 


HB (red clump) 


split 


broadened 



the value of the point, and with the sigma of the point, and we 
calculated the slope A for each of a million random samplings 
from these Gaussians. The million values of A had a distri- 
bution that was closely Gaussian; Table 5 gives the resulting 
slopes (per arcmin) and their sigmas. In the last column of the 
table is the ratio of \A\ to cr, which is an index of how reliably a 
slope diff'ers from zero; a 2cr result is of course only marginal, 
while a 3cr result is quite significant. Thus the radial gradient 
in Fig. 19 is marginal, that in Fig. 20 dubious, and those in 
Figs. 21 and 22 highly significant. 

5. Discussion 

5.1. Summary of the observational evidence of multiple 
stellar populations 

We have analyzed a large sample of proprietary (GO- 1 1739) 
and archival images of the two massive, metal-rich Galactic- 
bulge GCs NGC 6388 and NGC 6441, have used proper mo- 
tions to remove most of the field stars, have corrected all our 
CMDs for differential reddening, and have made a thorough 
multicolor analysis of the stellar populations of the two clus- 
ters. 

In both clusters we found significant evidence for the pres- 
ence of at least two stellar generations, which manifest them- 
selves as broadened or multiple sequences in all the main evo- 
lutionary branches of the CMDs. Table 6 collects the main 



26 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 O.B 1 

™F606W -'^^^FBUW ™ F606W FB14VI 



Fig. 23. — The horizontal-branches of NGC 6388 and NGC 6441 in the mpeoew vs. mpeoew - '«f814w CMD. We identified by hand four 
regions, colored in red for the RC, and green, azure, and blue for three distinct segments of the blue HB, from brightest to faintest. 
Internal labels give the number of stars in each of these regions, and the fraction of the total HB population that it constitutes. 



results on the multiplicity of the various sequences. These re- 
sults are distilled from a quite complex data set, whose most 
important aspects are summarized next. 

The MS of NGC 6388 is analyzed with various filter com- 
binations in Figs. 9-11. Most of the CMDs show a broadening 
of the MS, but none of them succeeds in splitting it into dis- 
tinct sequences as instead happens for NGC 6441 when using 
the ;7Zf39ow - '«f814w color, as shown in Figures 6(g) and 21. 
The separation of the two MSs of this latter cluster as a func- 
tion of color baseline is illustrated in Fig. [12] Figure |8] further 
confirms the reality of these results for both clusters, showing 
that neither of them could result from differential reddening 
(for which all CMDs have been corrected). Note, moreover, 
that the average reddening of NGC 6388 is lower than that of 
NGC 6441 (see subsection 2.5), hence it should have been eas- 
ier to recognize two MSs (if present) compared to the case of 
NGC 6441. We conclude that the MS structures of the two 
clusters are different. 

The SGBs of the two clusters are examined in Fig.[T3]which 
shows a clear split for the SGB of NGC 6388 in panel (a) but at 
best a doubtful one for NGC 6441 in panels (c) and (d), where 
the SGB components apparently overlap; Figures 14 and 15, 
respectively, explore the SGB for each cluster in multicolor 
domains. Again, we note this more subtle difference between 
the two clusters, with two distinct SGBs in NGC 6388 and just 
a broadened SGB in NGC 6441 (the opposite situation com- 
pared to the MS). 

The structure of the RGB of the two clusters is shown in 
Fig. 16. In both cases there is a hint of a double RGB, more 
evident in the upper part of the RGB of NGC 6388 and in the 
lower part for NGC 6441. The evolutionary connection be- 
tween the MS and the RGB is somewhat ambiguous in the case 
of NGC 6441: the population ratio seems to indicate that bMS 



and fRGB stars belong to the same population, but the radial 
gradient (which complicates the calculation of the population 
ratios) suggests instead to connect the bMS to the bRGB. For 
NGC 6388 it appears that bSGB and bRGB stars belong to the 
same population. 

For the HB we have thus far examined only the red clump; 
we take up the extended blue part of the HB in the following. 
The split of the RC of NGC 6388 can be seen in the two-color 
diagrams in panels (a) and (c) of Fig. 17, while for NGC 6441 
the broadening is seen in most of the panels of Fig. 18. There- 
fore there is some evidence that the RC is also a mixture of the 
first and second generations. However, what makes these two 
clusters unique, among the 150 globular clusters of the Milky 
Way, is their blue HB extension in spite of their high metal- 
licity. Furthermore, their blue HBs are remarkably similar, as 
illustrated in Fig. |23] Note also the similarity of the relative 
number of stars in the distinct sections into which we have di- 
vided the HB. This similarity is particularly surprising given 
that — as is well known — the distribution of stars on the HB 
is extremely sensitive to many different parameters (age, com- 
position, RGB mass loss). Thus, the differences in the structure 
of the MS, SGB and RGB of the two clusters demonstrate that 
there must be a sizable difference in some of those parameters, 
and yet the HBs turn out to be very similar. This conundrum is 
further discussed in the last subsection. 

5.2. The helium abundance of the secondary populations 

Photometric differences in the CMD between multiple stel- 
lar populations in GCs have been attributed to differences in 
helium (e.g., Bedin et al. 2004; Norris 2004; D'Antona et al. 
2005; Piotto et al. 2005, 2007) or in light-element abundances 
(e.g., Sbordone et al. 2011). Different photometric bands can 
enhance different effects of the chemical pattern of the popu- 
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Fig. 24. — Points from panel (d) of Fig. [12] on which we have 
superimposed values of Acolor that follow from synthetic col- 
ors calculated with our best-fit value of the helium abundance 
ofthebMS ofNGC 6441. 

lation. As noted in the introduction, the presence of helium- 
enriched stars in these two clusters was first proposed based 
on the exceptional morphology of their HBs and on the prop- 
erties of their RR Lyrae variables. Now, the discovery of a 
double main sequence in NGC 6441 lends support to this inter- 
pretation, and favors a scenario in which the helium-rich stars 
belong to a second stellar generation, as opposed to the possi- 
bility that they would owe high helium to deep-mixing on the 
upper RGB, as proposed by Sweigart & Catelan (1998). 

Our studies of 47 Tuc (Milone et al. 2012b), NGC 6397 
(Milone et al. 2012c), and NGC 6752 (Milone et al. 2010) 
demonstrate that light-element variations should not signifi- 
cantly affect the »Jf606w — '«f814w or wipsuw — wJfi60w col- 
ors of MS stars, so we explore the possibility of helium be- 
ing the chief cause of the observed MS split in NGC 6441 (or 
of MS broadening, depending on the specific bands). To this 
end, we calculated synthetic colors for the two MSs, adopt- 
ing a primordial helium 7=0.26 for the rMS, while trying for 
the bMS Y values between 0.27 and 0.36, in steps of 0.01. 
We assumed for each sub-population the same iron abundance 
([Fe/H]=-0.55 Harris 1996, 2010 edition, see also Carretta et 
al. 2009b, Origlia et al. 2008). Note that a difference of 0.1 
dex in [Fe/H] (which represents the present uncertainty in the 
metallicity difference between the two clusters) in the calcu- 
lated models does not change our conclusions. 

In previous papers we have used a similar approach to es- 
timate helium differences between GC sub-populations, but 
there we knew the light-element abundance of the two MSs 
and here we do not. Hence, we assume the abundance of C, N 
and O to be the same, in spite of likely differences in at least 
their relative proportions between the first and second gener- 
ations. In practice, we assumed [C/Fe]=-0.45, [O/Fe]=0.27 
(Origlia et al. 2008), and solar N abundance. In both cases we 
used [afFe]-Q.3 (Origlia et al. 2008). Neglecting possible dif- 
ferences of CNO proportions should not affect our conclusion 
about helium because, as noted above, the F606W, F8 14W and 
F160W filter bands are almost completely insensitive to C, N, 



and O, while only F390W is expected to be somewhat affected 
by abundance variations in these light-elements. 

We adopted the synthetic colors from the BaSTI isochrones 
(Pietrinfemi et al. 2004) in order to compare observations and 
synthetic colors, and determined Teff and log^ for MS stars 
at mF8i4w=21.5. The ATLAS12 code (Kurucz 2005; Castelli 
2005; Sbordone et al. 2007) was then used, which allows us 
to work with arbitrary chemical compositions. We performed 
spectral syntheses from 3,400 A to 20,000 A by using the SYN- 
THE code (Kurucz 2005), and integrated the synthetic spectra 
over the transmission curves of the four filters. From the re- 
sulting fluxes we formed, for each value of Y, the differences 
Acolor = mx — mpguw- 

Our result is shown in Fig. |24] We found that the value 
ibMS=0.32 optimizes the agreement between models and ob- 
servations. The overall good agreement suggests that helium 
can indeed be the main cause of the observed split, with bMS 
stars being He-enhanced by Ay~0.06 relative to rMS stars. 
Best-fit model and data are less in agreement for the F390W 
point, as expected because this band-pass includes CN and CH 
molecular bands, hence the flux through this band is sensitive 
to CNO variations (e.g., Milone et al. 2012b). 

In the case of NGC 6388 there is no evidence of a split in its 
MS, which is clearly broader than expected from photometric 
errors, but the color spread (in ;7Zf39ow -"Zf814w) is smaller than 
that of NGC 6441. However, as previously emphasized, the 
'Wf606w-'«f814w color baseline is more sensitive to helium vari- 
ation, but is insensitive to variations of CNO. In Figure |25] we 
compare the mF606w vs. otf606w - '«f814w CMDs of NGC 6388 
and NGC 6441. Panels (a) and (d) are the CMDs themselves 
(from the overlapping regions of GO-11739 and GO-10775, 
proper-motion-selected members only for NGC 6441). Panels 
(b) and (e) show the rectified MSs, extending ~3 magnitudes 
below the TO. The rectified Acolor distribution, in five magni- 
tude intervals, is shown in panels (c) and (f). 

The color distributions appear to be symmetric in NGC 6388 
but slightly skewed bluewards in NGC 6441. Table 7 gives 
the MS breadths, from fitting Gaussians to the distribution of 
Acolor The first column identifies the panel in the figure, 
and the final column gives the intrinsic breadths, calculated as 
the quadrature difference between the preceding two columns. 
These breadths are also plotted in Fig. |26]as a function of the 
'«F606w - "^6()6w magnitud^. We conclude that the intrinsic 
WJF606W - WZF814W color dispersions of the two MSs are quite 
similar, so that despite the different manifestation of the two 
populations, at least their helium distributions must likewise 
be quite similar. 

Though a complete and thorough analysis and interpretation 
of the numerous CMDs presented here is deferred to a future 
paper (Cassisi et al. in preparation), we wish to present here 
some further remarks on the estimates of the helium abundance 
in the stellar populations of NGC 6388 and NGC 6441 . 

The F~0.32 that we found for the bMS stars in NGC 6441 
from the MS separations in the various color baselines can 



To choose the TO value in hif606W we derived a fiducial fine along the MS and 
took it's bluest point. 
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Fig. 25. — Illustration of the procedures used to determine the MS intrinsic breadth of the two clusters in /iipgoew ~ w^fsmw- Panels (a), 
(b), and (c) refer to NGC 6388, Panels (d), (e), and (f) to NGC 6441. See text for details. 



be compared with the He content estimated by Caloi and 
D'Antona (2007) and by Busso et al. (2007) from their analy- 
sis of the HB morphology and of the pulsational properties of 
the RR Lyrae stars in this cluster Busso et al. were able to 
reproduce the blue and tilted HB morphology of NGC 6441 by 
assuming a He abundance yaiO.35 for the second-generation 
stars, only ~0.03 higher than our estimate from the MS split. 
At the same time Caloi & D'Antona suggested that a fraction 
of the NGC 6441 stars would have a He content as high as 
7=0.38-0.40 in order to account for the extreme blue part of 
the HB. Part of this discrepancy with our present result and 
that of Busso et al. could be due to the use of optical bands by 
Caloi & D'Antona, which surely are not optimal for studying 
the very hot HB extension (see Dalessandro et al. 2011). More- 
over, in order to match synthetic to observed HBs, both Busso 
et al. and Caloi & D'Antona appeal to a substantial spread of 
helium from 7^0.25 in the first generation up to the their max- 
imum estimated values (0.35 and 0.40, respectively). Our re- 
sult, with two distinct MSs in NGC 6441 suggests instead that 
(at least in this cluster) there are basically two distinct values 
of the heUum abundance, with a clear gap in between. 

Furthermore, there is yet another puzzling inconsistency be- 
tween the results of those two investigations based on the HB 
and the results of the present paper According to Caloi & 
D'Antona only -38% of the stars in NGC 6441 would have 
primordial He content {Y ~0.25), whereas according to Busso 
et al. only 15-16% of the stars should belong to the second 
generation. Though in different ways, both results are at vari- 
ance with what we obtained from the two MSs of this cluster, 
i.e., the majority of stars (~ 65%) are on the red (presumably 
first-generation) MS and only ~ 1/3 appear to belong to the 
second generation. We shall come back to this issue in the next 
subsection. 

At least part of these discrepancies can be traced to the as- 
sumption made by Busso et al. that the whole red clump con- 
sists of first-generation stars, with second-generation stars con- 
fined to the blue HB extension; indeed, the RC represents ~ 12- 
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Fig. 26. — MS breadth in ffipeoew - wJfsmw as a function of the 
magnitude mpeoew below the TO. Blue for NGC 6388, red for 
NGC 6441. 
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Intrinsic MS Breadth in mpeoew • 
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15% of the total HB population. We have found instead that the 
RC is a mixture of first- and second-generation stars. It is then 
possible that all second-generation stars arrive on the HB at 
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the RC, spend there roughly 50% of the HB lifetime, and then 
evolve into an extended blue loop for the remaining ~50% of 
the time. In fact, as emphasized by Caloi & D'Antona, the 
evolution of the helium-rich and metal-rich HB models may be 
confined to the RC region, or, after spending some time on the 
RC, the stars may experience an extended excursion reaching 
very blue colors for differences in total mass of order of just 
~ 0.01 Mq (Sweigart & Gross 1978). For this reason, it is difli- 
cult to precisely pinpoint the helium abundance from the mere 
extension of the HB. 

Finally, we note that the MS breadth in the mpeoew - ffiF8i4w 
color increases with decreasing luminosity (see Fig. 26 and Ta- 
ble 7) in a fashion that apparently cannot be accounted for by 
the increasing photometric errors. Here, we simply observe 
that an increase of the MS spread in the 2-3 magnitude inter- 
val below the TO is also observed in the MSs of a> Cen (King 
et al. 2012) and NGC 2808 (Piotto et al. 2007). We shall return 
on this observational evidence in a future paper. 

5.3. An NGC 6388-NGC 6441 conundrum? 

As already emphasized, these two clusters are unique in 
having a major blue horizontal branch extension in spite of 
their high metallicity and in the properties of their RR Lyrae 
variables. Furthermore, their HBs are remarkably similar in 
shape and extension and even in the relative star numbers along 
the branches (see Fig.l23ll. The two clusters also have in com- 
mon the presence of two population components whose dis- 
tinctness is evident in various parts of their CMDs, with the 
population that is (slightly) more concentrated to the center 
likely being the second generation. The outstanding anomaly, 
however, is that whereas the MS of NGC 6441 is clearly split 
into two branches (in mFsgow-wiFSMw color), that of NGC 6388 
is single, though broadened. Given the extremely high sen- 
sitivity of the HB morphology to the helium abundance, we 
conclude that the helium difference between the two stellar 
generations must be similar, in the two clusters. The results 
presented in Fig.|25]and|26] showing the close similarity of the 
MS of the two clusters in the mF606w -'Wf814w color, strengthen 
this conclusion further. How can these clusters exhibit almost 
identical HBs while having very different MS properties in 
WF390W - '«F8i4w color? In other words, if the second gener- 
ation populations in the two clusters have nearly identical he- 
lium abundances, why is the MS split in NGC 6441 but not in 
NGC 6388? Clearly, something else must differ in order to re- 
duce below detection the MS split of NGC 6388, or to amplify 
that of NGC 6441 . This is the conundrum of what is otherwise 
a pair of twin clusters. 

One possibiUty is metalUcity [Fe/H], such that in NGC 6388 
the helium-rich population is also somewhat metal-enriched, 
while it is not in NGC 6441. This possibility appears to be 
excluded by the high-resolution spectroscopy of Carretta et al. 
(2007), according to which RGB stars in NGC 6388 have the 
same [Fe/H] within ~0.01 dex. Their sample consists of just 
seven stars, however, so that the chance of their having missed 
the minority population is not entirely negligible. However, 
among these seven stars there appears to be a spread in oxy- 
gen, sodium, and aluminum abundances, which should indi- 



cate that, though small, the sample actually includes stars from 
both the first and the second generation. This is confirmed also 
by a lower-resolution study of a larger sample of stars (Car- 
retta et al. 2009a). The option that the seven stars of Carretta 
et al. all belonged to a single population does not appear to be 
supported by observations. 

As emphasized above, the main difference between the two 
clusters arises when the F390W filter comes into play, and 
within its passband strong CN and CH bands are included (see 
Fig. 32 in Milone et al. 2012b). A more attractive option is 
therefore that the secondary population in NGC 6388 is CNO- 
enhanced compared with the first generation, though this re- 
quires a chemical enrichment scenario that is quite ad hoc. In 
this way excess CN and/or CH molecular blanketing in the UV 
would make redder the MS of the secondary population in this 
cluster. A similar effect could also be produced by the C:N;0 
proportions in the secondary population of NGC 6388 being 
different from those in the secondary population of NGC 6441 . 

We defer to a future paper a more quantitative exploration 
of these alternatives (Cassisi et al., in preparation), but we note 
that further progress will require an extensive spectroscopic ex- 
ploration of these clusters at high spectral resolution, in partic- 
ular having photometrically identified the targets as belonging 
to one or the other sub-population. More extensive HST obser- 
vations, especially at ultraviolet wavelengths, would also help 
in solving this conundrum. 
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